Effect of matrix parameters on mesoporous matrix based quantum computation 
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We present a solid state implementation of quantum computation, which improves previously proposed op- 
tically driven schemes. Our proposal is based on vertical arrays of quantum dots embedded in a mesoporous 
material which can be fabricated with present technology. We study the feasibility of performing quantum 
computation with different mesoporous matrices. We analyse which matrix materials ensure that each individ- 
ual stack of quantum dots can be considered isolated from the rest of the ensemble-a key requirement of our 
scheme. This requirement is satisfied for all matrix materials for feasible structure parameters and GaN/AIN 
based quantum dots. We also show that one dimensional ensembles substantially improve performances, even 
of CdSe/CdS based quantum dots. 



I. INTRODUCTION 

Solid state based hardware is considered a promising can- 
didate for implementing quantum computing 1 . This is due to 
the potential scalability of the involved nanostructures, and 
to the relative ease of integrating them with more traditional 
electronic components, e.g. for input/output functions. Some 
interesting proposals rely on spin and exciton qubits con- 
fined in semiconductor self-assembled quantum dots (QD's) 2 , 
which can be manipulated using ultra fast laser pulses 3,4,5 . 
In these schemes qubit coupling, which is necessary to per- 
form universal quantum computation, can be realised through 
exciton-exciton direct Coulomb interaction 6,7,8 . The presence 
of an exciton in a QD produces a biexcitonic shift in the 
ground state excitonic energy of a nearby QD. By driving a 
qubit at this shifted frequency, conditional operations can be 
performed 6,7 . 

Currently two major problems for QD based quantum com- 
putation are addressability of the single qubit and size poly- 
dispersity in self-assembled QD ensembles. In order to carry 
out practical quantum computation it is necessary to be able 
to address individual qubits. This poses a problem for op- 
tically driving the response of self assembled quantum dots 
(such as the ones grown by Stranski-Krastanow techniques): 
in these dot ensembles the size of each dot is one order of mag- 
nitude smaller than that of the laser spot addressing it. The 
use of energy selective methods on isolated stacks of quan- 
tum dots (quantum registers) has been proposed to solve this 
problem 6,7 . It is still experimentally difficult to control QD 
size and position in a satisfactory way, and QD vertical stacks 
tend to form at random positions. Additionally the size of the 
QDs within the stacks is mainly due to strain propagation and 
it is hardly controllable, resulting in the practical difficulty of 
creating the desired sequences of energy selectable excitonic 
transitions. To overcome these difficulties we recently pro- 



posed to use mesoporous matrices to create regular arrays of 
QD qubits suitable for quantum computation 9 . The redundant 
encoding typical of the chosen hardware protects the compu- 
tation against gate errors and the effects of measurement in- 
duced noise. We have discussed that quantum computation 
can be efficiently implemented using this hardware for both 
III-V and II- VI materials in a TiC>2 matrix. In this paper we 
analyse how the matrix material affects the implementation of 
quantum computation. 



II. COMPUTATIONAL HARDWARE 



We consider alternating layers of two semiconductors with 
widely different band gaps deposited in a mesoporous matrix. 
This provides a stack of QDs (qubits), sandwiched between 
the larger band gap material (barriers). The resulting system 
is depicted in Fig. 1(a) and consists of an array of identical, 
hexagonally-packed stacks of quantum dots (columns). The 
band structure within each column is sketched in Fig. 1(b). 

We consider semiconductors which assume wurtzite crystal 
structure for both the barriers and qubits. This leads to strong 
intrinsic electric fields 10 within each column, (see Fig. 1(b)). 
The built-in electric field enhances coupling between excitons 
in neighbouring quantum dots within a stackL. Quantum com- 
putation can then be carried out by using sequences of laser 
pulses, as described in Ref. 3 . Under the influence of the same 
laser pulse, each column within the laser spot will act as an 
independent replica of the same computational array. This in- 
trinsic large redundance protects the computation even in the 
presence of errors and measurement induced noise 9 . 
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FIG. 1 : (a) The proposed system, an ensemble of stacks of alternating 
QDs and barriers embedded in a mesoporous matrix, (b) The band 
structure of each individual stack including the intrinsic field of the 
materials. 



III. THEORETICAL MODEL 

Each individual stack is modelled as a column of cylindri- 
cal quantum dots with the same radius. To calculate the biex- 
citonic shifts, the confining potentials are approximated as 
parabolic potentials with the same characteristic widths as the 
expectation values V < z 2 > and V< r 2 > of the cylindrical 
dots in the stack. Nearby quantum dots (qubits) are coupled 
by the biexcitonic shift AE between ground state excitons. 
In order to calculate interaction energies between qubits we 
approximate the biexcitonic shift AE as 



\AE\ 



(^1^2 | U C |^1 ^2) 



(1) 



where ipi(r ie ,r ih ) = ^ie{r ie )^i h {r ih ) is the wavefunction of 
the exciton in QD^ in the single particle approximation, and 
Uc is the Coulomb interaction between the two excitons. In 
order for the scheme to be feasible two conditions must be 
satisfied. The biexcitonic shift between neighbouring qubits 
within a stack must be large enough, so quantum computation 
can be performed within the stack on time scales much shorter 
than the relevant decoherence times, and the inter- stack inter- 
actions must be negligable so that each QD array can be con- 
sidered isolated. 

The biexcitonic shift between excitons within a stack is ex- 
ploited to perform two-qubit gates using multicolour trains 
of laser pulses 6 ^. For performing operations on picosecond 
time scales - which is essential due to the relatively short ex- 
citonic decoherence times - AE must be of the order of a few 
meV. When choosing the correct parameter range however, 
additional factors must be taken into consideration 6 . Larger 
AE can be induced by increasing the height of the quan- 
tum dots. This increases the excitonic dipole moments un- 
der the influence of the intrinsic electric field but reduces the 
oscillator strength, i.e. the efficiency of driving the compu- 
tation through coupling to laser pulses. Finally the barrier 
width must be large enough to ensure that single particle tun- 
nelling between quantum dots is negligible on the relevant 
time scales. 

If all the stacks in the ensemble were to compute correctly 
-i.e. no computational errors - inter- stack interactions would 



only renormalise the excitonic energies by the same amount 
for each stack, which could be easily accounted for in the 
computational scheme. However the event of computational 
failure in a certain stack would induce a local, unwanted, shift 
of the exciton energies in neighbouring stacks. Each individ- 
ual stack will interact with a certain number of such 'faulty 
computing' stacks. Therefore the magnitude of this unwanted 
shift is the sum of the interaction energies of a resident exciton 
with an exciton in each 'faulty computing' stack |A£^ er |. 
The inter-stack interaction between qubits is negligible, if it is 
much smaller than the interaction between excitons within a 
column, i.e. |A£^ er | << \AE\. The stacks are arranged in 
a hexagonal structure, which can be represented as a series of 
concentric hexagonal shells surrounding each stack. The i th 
shell consists of 6i stacks. The number of expected computa- 
tional failures in the i th shell will therefore be 6i • (1 — p), 
where p is the probability of a successful computation for any 
individual stack. We estimate the total energy shift due to 
interactions with all the 'faulty computing' stacks in the en- 
semble as: 



ensemble 



\AEt 



^ 6i ( 1 - p) AEi n ter ) 



(2) 



where = (r max + r m i n )/2 is the average radius of the i th 
hexagonal shell, r max and r m ^ n being the maximum and min- 
imum distance of the shell from the central stack, and AE inter 
is calculated according to Eq. (Q]) considering dots in different 
stacks. 



IV. INFLUENCE OF MATERIAL PARAMETERS AND 
DIMENSIONALITY 

We have shown 9 that for typical III-V and II- VI materi- 
als embedded in a TiC>2 matrix, the computational hardware 
can be designed such that the biexcitonic shift between qubits 
within a stack is large enough to perform quantum computa- 
tion, and that the inter- stack interactions are small enough for 
the hardware to be considered an ensemble of isolated stacks. 
For a GaN QD and a A1N barrier, or for a CdSe QD and a CdS 
barrier, each of radius 5nm, there exists a range of QD heights 
and barrier widths which lead to a sufficient biexcitonic shift 
between neighbouring qubits (AE > 3me V), satisfactory os- 
cillator strength (/i > and a sufficiently long 
tunnelling time between qubits (r > Ins). The tunnelling 
time is estimated by taking the inverse of the tunnelling rate 
as expressed in Ref J^. 

The interaction between neighbouring stacks can be de- 
creased by increasing the thickness of the matrix walls (see 
Fig. 2.), which can be made up to « 6 nm thick before the 
order of the array is lost. The magnitude of the interaction 
with faulty computing stacks depends also on the probability 
p of each stack computing successfully (see Eq (2)). It was 
found 9 that for the GaN/AIN system with TiC>2 matrix, and 
matrix wall thicknesses of 6nm, each stack could be consid- 
ered isolated from the rest of the ensemble for all p. For a 
CdSe/CdS hardware it was found that for matrix walls 6 nm 
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FIG. 2: Interstack qubits separated by the matrix wall. 
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FIG. 3: R versus the number of shells included in the calulation, for 
the GaN/AIN system in a TiCb matrix with 2 nm walls, and p — 0.83 
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FIG. 4: R(oo) for a GaN/AIN system versus the dielectric constant 
of the matrix for different p. a) matrix walls 2 nm thick and b) matrix 
walls 6 nm thick matrix walls. 
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thick the stacks could be considered isolated for p ~ 0.86, 
and similarly for a wall thickness of 4.5 nm for p ~ 0.91. 

Figure 3 shows the ratio R = \AEj°l er \/\AE\ against the 
number of shells of neighbours included in the calculation, for 
the GaN/AIN system in a Ti0 2 matrix (e T io 2 ~ 100) 12 with 
2 nm walls, and p = 0.83. R(oo) is the number at which this 
ratio saturates. For the interstack interaction to be negligable it 
must be an order of magnitude smaller than AE i.e. R(oo) < 
0.1. For the case in Fig. 3. R(oo) « 0.066, therefore the 
interstack interactions can be neglected and each stack in the 
ensemble can be considered isolated. 

Ti02 was initially proposed as the matrix material due to 
its high dielectric constant, resulting in high screening be- 
tween stacks. However experimentally other materials may 
be used to produce mesoporous matrices. In the following 
we will analyse how different materials, characterised by dif- 
ferent dielectric constants, would affect the feasibility of our 
scheme by increasing or reducing the inter- stack interaction. 
Fig. 4. shows R(oo) for a GaN/AIN system against the di- 
electric constant of the matrix for different p for matrix walls 
a) 2 nm thick and b) 6 nm thick. It can be seen that the con- 
dition R(oo) < 0.1 still holds for materials with dielectric 
constants much lower than T1O2 . For example a much more 
standard material for the matrix is S1O2 which has a dielec- 
tric constant e w 4. Fig. 4. shows that for a S1O2 matrix with 
walls as thin as 2 nm andp ~ 0.9 the inter- stack interaction is 
already negligable. For 6 nm walls S1O2 is a suitable matrix 
for p ~ 0.83. We underline that p = 0.83 corresponds to, on 
average, one nearest neighbour failing. The situation is differ- 
ent for a CdSe/CdS based system. As shown in Fig. 5, when 



FIG. 5: R(oo) for a CdSe/CdS system versus the dielectric constant 
of the matrix for p = 0.83, for a) matrix walls 6 nm thick and b) 
3 nm thick matrix walls. 



p = 0.83, materials with dielectric constant e > 130 would be 
needed to achieve an ensemble of isolated stacks. 

One possible solution to decrease the interaction between 
stacks for the II- VI system would be to reduce the dimension- 
ality of the ensemble array. By reducing the 2D array of stacks 
to a ID chain the interaction energy between a single stack and 
the 'faulty computing' stacks in the ensemble is reduced to 



\AE? n t ter \ = 2(l-p) 



ensemble 

E 

i=l 



AE inter (n). (3) 



Figure 6 shows the ratio R, against the number of neighbours 
included in the calculation for a T1O2 matrix and a chain of 
CdSe/CdS-based stacks. It is shown that for a wall thickness 
of 3 nm the chain can be considered isolated for p as low as 
0.76 (and for even lower p when the walls are thicker). This is 
a large improvement in respect to the results for the 2D array 
shown in Fig. 5b. Reducing the dimensionality of the array 
has not been experimentally achieved. If the conditions during 
manufacture could be changed to produce a rectangular lattice 
instead of hexagonal (this change alone would reduce the in- 
teraction energy between stacks), an effective sytem of ID 
chains could be realised by having one lattice constant much 
greater than the other. Such a configuration in fact could be 
seen as a group of independent chains of stacks. 



4 



0.25 




10 20 30 40 50 60 70 80 90 100 
i 

FIG. 6: R for a chain of CdSe/CdS stacks, against the number of 
pairs of nearest neighbours included in the calculation for p = 0.83, 
and for 3 nm thick matrix walls. 



V. CONCLUSIONS 

Our proposal is based on a mesoporous matrix which pro- 
vides an uncorrelated ensemble of computational arrays. This 
intrinsic redundancy improves the response of the system, 
when used as quantum computation hardware 9 . We have dis- 



cussed how the use of different matrix materials would in- 
fluence this property, namely by increasing/decreasing inter- 
stack interactions. It has been shown that, for GaN-based 
qubits, the stacks in the ensemble can be considered isolated 
for any reasonable matrix material, provided that the matrix 
walls are 6 nm thick and the probability of each isolated stack 
in the ensemble computing successfully is greater than a 0.83. 
The stacks in the ensemble can still be considered isolated for 
matrices with thinner walls as long as there is a higher prob- 
ability that each individual stack will compute correctly. It is 
noteworthy mentioning that, when p ~ 0.9, even a material 
with a relatively low dielectric constant such as Si02 could be 
used, and with matrix walls as thin as 2 nm. 

Finally we have shown that reducing the dimensionality of 
the ensemble also improves the response of the system, by re- 
ducing the influence of interactions between the members of 
the ensemble. This is particularly important for CdSe/CdS- 
based systems. In this case, for a T1O2 matrix with hexago- 
nally packed pores, to be able to regard the QD stacks as inde- 
pendent, the matrix walls must be at least 4.5 nm thick and the 
probability of individual success as high as p — 0.9 9 . How- 
ever the use of matrices with pores aligned along ID chains 
could allow the probability of success for isolated stacks to be 
as low as 0.76, with matrix walls as thin as 3 nm. 



1 For a general review on quantum computation, see e.g M. A. 
Nielsen and I. L. Chuang "Quantum Computation and Quantum 
Information" (Cambridge University Press 2000) 

2 P. M. Petroff, A. Lorke, and A. Imamoglu, Physics Today 54, issue 
5, 46 (2001). 

3 see e.g. "All-Optical Schemes for Quantum Information Process- 
ing with Semiconductor Macroatoms" I. D' Amico, E. Pazy, P. Za- 
nardi, E. Biolatti, R.C. Iotti, F. Rossi, F. Troiani, U. Hohenester 
and E. Molinari, in 'Semiconductor Macroatoms: Basic Physics 
and Quantum-Device Applications' (2005) (edited by Imperial 
College Press/World Scientific Publishing). 

4 A. Nazir, B. W. Lovett, S. D. Barrett, J. H. Reina, and G. A. 
Briggs, Phys. Rev. B 71, 045334 (2005). 

5 X. Li, Y. Wu, D. Steel, D. Gammon, T. H. Stievater, D. S. Katzer, 
D. Park, C. Piermarocchi, L. J. Sham, Science 301, 809 (2003). 



6 E. Biolatti, I. D' Amico, P. Zanardi and F. Rossi , Phys. Rev.B 65, 
75306 (2002). 

7 S. De Rinaldis, I. D'Amico, E. Biolatti, R. Rinaldi, R. Cingolani, 
and F. Rossi, Phys. Rev. B 65, 081309(R) (2002). 

8 E. Pazy, E. Biolatti, T. Calarco, I. D'Amico, P. Zanardi, F. Rossi, 
P. Zoller, Europhys. Lett. 62, 175 (2003). 

9 T. E. Hodgson, M. F. Bertino, N. Leventis, and I. D'Amico, J. 
Appl. Phys. 101, 114319 (2007). 

10 R. Cingolani, A. Botchkarev, H. Tang, H. Morkoc, G. Traetta, G. 
Coli, M. Lomascolo, A. Di Carlo, F. Delia Sala, and P. Lugli, 
Phys. Rev. B 61, 271 1 (2000). 

11 E. Pazy, I. D'Amico, P. Zanardi, and F. Rossi, Phys. Rev. B 64, 
195320 (2001). 

12 R. D. Shannon, J. Appl. Phys. 73, 348 (1993). 



